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We calculate the relative amplitude of orbital angular mo- 
mentum (0AM) entangled photon pairs from the spontaneous 
parametric down conversion. The results show that the am- 
plitude depends on both the two Laguerre indices I, p. We 
also discuss the influences of the mostly used holograms and 
mono-mode fibers for mode analyzation. We conclude that 
only a few dimensions can be explored from the infinite OAM 
modes of the down-converted photon pairs. 

PACS number(s): 03.67.Mn, 03.65.Ud, 42.50.Dv 



I. INTRODUCTION 

Quantum entanglement is a very important property 
of quantum mechanics. It is the foundation of quantum 
teleportation [1], quantum computation [2-4], quantum 
cryptography [5], superdense coding [6], etc. Up to now 
most of the theoretical discussions and experiments are 
focused on quantum states belonging to two-dimensional 
states, or qubits [7-10]. In recent years, the interest in 
multi-dimensional states, or qudits, is steadily growing 
for its promise to realize new types of quantum commu- 
nication protocols [11-13], and its properties in quan- 
tum cryptography better than qubits [12-15]. These 
theoretical disscussions [16-20] and experiments [21-25] 
about multi-dimensional states are mostly based on the 
orbital angular momentum (OAM) of the photons. It 
has been shown that paraxial Laguerre-Gaussian(LG) 
laser beams carry a well-defined orbital angular mo- 
mentum [26], and that the LG modes form a complete 
Hilbert set. The down-converted photons from sponta- 
neous parametric down-conversion (SPDC) are entangled 
in not only polarization, or spin angular momentum, but 
also OAM [22,24]. This provides a promise to explore 
multi-dimensional quantum state in one photon [22,18]. 

In this paper, we calculate the relative amplitude of 
OAM of the down-converted photons from SPDC and 
analyze the possible joint detection probability under 
the influence of the experiment elements. Our results 
show that the relative amplitude decreases almost expo- 
nentially with the growing of OAM. And only a few di- 
mensions can be explored from the infinite OAM modes 
of the down-converted photon pairs. In Section 2, we 



briefly introduce the LG mode, and calculate the rela- 
tive amplitude of every LG mode of the down converted 
photons from SPDC in detail. In theory, the relative am- 
plitude determines the joint detection probability. How- 
ever, in practical experiments, the computer generated 
holograms and the mono-mode fibres will inevitably in- 
fluence the joint detection probability. We discuss the 
mode analysis after the computer generated hologram in 
Section 3. Section 4 analyzes the detection efficiency of 
mono-mode fibre for every LG mode. Section 5 presents 
the possible joint detection probability, when the effects 
of the holograms and mono-mode fibres are both in- 
cluded. The last section is the conclusion. 



II. SPONTANEOUS PARAMETRIC DOWN 
CONVERSION AND OAM 

It is well known that photons can carry both spin angu- 
lar momentum and OAM [26]. Spin angular momentum 
is associated with polarization and OAM with the az- 
imuthal phase of the electric field. The normalized LG 
mode is given in cylindrical coordinates by 
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where Lp(x) are the associated Laguerre polynomials. 
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and the standard definitions for Gaussian beam parame- 
ters are used: 

uj{z) — ujQy^l + (z/zfl,)2 : spot size, 

R{z) — z{\ -\- [zrI zY) : radius of wavefront curvature, 

tp{z) = arctan(z/z/j) : Gouy phase, 

zr = ^kujQ : Rayleigh range. 

Wo is the beam width at the beam waist, the index I 
is referred to as the winding number, and (p + 1) is the 
number of radial nodes. If the mode function is a pure 
LG mode with winding number I , then every photon of 
this beam carries an OAM of lU. This corresponds to an 
eigenstate of the OAM operator with eigenvalue Ifi [26]. 
If the mode function is not a pure LG mode, the state is 
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a superposition state, with the weights dictated by the 
contribution of the ^th angular harmonics. 

At the beam waist {z = 0), the LG mode can be writ- 
ten as 
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In the following calculation, we use this equation be- 
cause in the experiment we always manipulate the light 
at its beam waist. 

In the SPDC process, a thin quadratic nonlinear crys- 
tal is illuminated by a laser pump beam propagating in 
the z direction, with wave number kp and waist loq. The 
generated two-photon quantum state is given by [20] 
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where {li,pi) corresponds to the mode of the signal beam 
and (^2,^2) the mode of the idler beam. The probability 
amplitude C^p\% is given as [16,17,19,20] 
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where r± is the radial coordinate in the transverse X~Y 
plane, $(rj_) is the spiral distribution of the pump beam 
at the input faced of the crystal, LGp{r±) is the spiral 
distribution of the LG mode beam at the same plane. 
The weights of the quantum superposition are given by 

^pVp2 ^ I^pi'p2l • ^^"^ ideal joint detection proba- 
bility for finding one photon in the signal mode 
and one photon in the idler mode {l2,P2)- 

Consider the case that the pump beam is in a pure LG 



mode iGpo with po = 0. 
can be written as 
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mode light at 2; = 
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Substitute LGq (p, (p) for ^{r±) into Eq. (5), and use the 
QAM conservation law in SPDC [16,22]: 
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where li.l2,lo are the winding numbers of signal beam, 
idler beam and pump beam respectively, we can achieve 
the probability amplitude 
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In the case = 0, or the input beam is Gaussian mode 
light, Eq. (8) can be simplified as(^ > 0) 
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It can be easily proved that C^'^.p^ = C'pi'ps = C'ps.pi = 
f-i.i 

P2,P1 ■ 

Table 1 gives the relative value for pi, p2 =0,1,2 
and I = 0,1,2. We can also illustrate the dependence of 
the relative probability amplitude Cp^j,^ on pi,p2, and 
I, with Fig. 1(; = and pi,p2 = 0,1,2,3,4) and Fig. 
2(pi = p2 = and ? = 0, 1, 2, 3, 4). 

From the above table and figures, we can see that the 
probability amplitude decreases very rapidly with the 
growing of pi, p2 and I. We then just consider the cases 
with pi, P2 = 0, 1, 2 and I = 0, 1, 2 when Po = Iq = 0. 
In papers [22,24], they also just consider the cases with 
I = 0,1,2. 

If additionally assume pi = P2 = as in the previous 
works [22,24,25,20], we can obtain: 
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This result can also be achieved from the Eq. (14) of the 
paper [20], when the condition li = —I2 = I is assumed. 

But till now, this assumption has not been proven ei- 
ther in theoretical discussions [16,20] or in experiments 
[22,24]. We will discuss the two cases separately in Sec- 
tion 5. Before proceed, we analyze the two main ex- 
periment elements, computer generated holograms and 
mono-mode fibre, which will unavoidably affect the de- 
tected relative probability amplitude. 



III. COMPUTER GENERATED HOLOGRAMS 
AND THE MODE ANALYSIS 

In most of the recent experiments [21-25], the authors 
always use computer generated holograms to transform 
Gaussian mode light into other LG modes light, or change 
the winding number of LG mode light. It is a kind of 
transmission holograms with the transmittance function: 

1 27r 
T{p,ip) = exp{i5— mod{l(fi — pcos</5, 2?:)), (11) 

where 6 is the amplitude of the phase modulation, A = |2L 
is the period of the grating at a large distance away from 
the fork, kx is the x component of the simplest reference 
beam's wave vector. Corresponding to the diffraction or- 
der m, the hologram can change the winding number of 
the input beam by Aim = 'ml- The diffraction efficiency 
depends on the phase modulation 5. When S = 27r, al- 
most 100% of the incident intensity is diffracted into the 
first-order. 
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However, even the input beam is a pure LG mode light, 
the output beam after the hologram is not a pure LG 
mode light. The output light will be the superposition 
of the various LG modes with the same I, and different p 
[27]. In addition, the beam waists of the input and out- 
put beam will affect the weights of different components 
of the output beam. Assume the input beam and the 
output beam have the same waists coq, thus the complex 
expansion coefficients of the decomposition of the mth 
diffraction order can be calculated as: 

api= j j {LG''^,^{p, if) exp{-im^r cos if))* 

T{p,ip)Ein{p,ip)pdpd<fi, (12) 

where Ein{p,(p) = iGp\(p, Consider the first-order 
diffraction, or m = 1, Eq. (12) can be then rewritten as 

/ j {LG\{p,^)r eM-i^l'f))LG';^{p,^)pdpd^, 

(13) 

where A/ = /i — is the winding number changed by 
the hologram. The relative weight of the output modes 
is given by 
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As the mono-mode fibres can only detect the photons 

with I — 0,we consider the case that the output light is in 

the modes with /'^ = 0. Then -P^^'p/ can be simplified as 



P^^ , , and P-^' , = P , 



In most of the experiments 

[22-25], only the holograms of = 1 or 2 are employed. 
Table 2 and Table 3 give the relative weight of different 
modes after the computer generated hologram with Al = 
1 and 2. 

Prom Table 2 and Table 3, wc can see most of the input 
mode LGp'^ is converted into LGp and LGp_^_i. Thereby, 
we only consider the case that the output light is in the 
modes with p' = 0, 1, 2, 3. 



IV. MONO-MODE FIBRE AND DETECTION 
EFFICIENCY OF THE OAM MODES 

It is known that only one mode of light can transmit 
in the mono-mode fibre: HEn mode. And in practical 
calculation, we always use Gaussian mode to replace the 
HEii mode. The Gaussian mode is 



E{p)=E{0)eM-pV^^), 



(15) 



where d = 2(jj is the Mode Field Diamcter(MFD) of the 
fibre, E{0) is amplitude of field at the fibre center. For 
LGp mode light, the detection efficiency is given as 
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Obviously if Z 7^ 0, then Qi p = 0. For the case I = 0, 
Eq. (16) can be simplified as: 
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To calculate the relative joint detection probability of 
the down-converted photons from SPDC, we only need 
the relative detection efficiency of the LG", LG\ , PG" 
and PGg, but not their absolute efficiency. Assume the 
waist size of the input beam is adjusted equal to d/2. 
Then when the detection area is miich more larger than 
the cross-section of the input light , only LGq mode light 
can be detected. But in practice, the detection area is de- 
termined by the fibre diameter. To simplify calculation, 
we further assume that the detection area is a round area 
with diameter equal to the MFD. Thus the integral for 
p is from to to. With these assumptions, the relative 
efficiencies for p = 0, 1, 2, 3 can be written as 

Qo ■■ Qi ■■ Q2 ■■ Qi = I ■■ 0.263 : : 0.036 (18) 



V. RELATIVE JOINT DETECTION 
PROBABILITIES OF OAM ENTANGLED 
PHOTONS FROM SPDC 

With the above discussions about the influence of com- 
puter generated holograms and mono-mode fibres, we 
now consider the joint detection probability for OAM en- 
tangled photons generated from an experimental set-up 
similar to the work [22] . The joint detection probability 
can be written as 

2 2 Pl+lp2+l 



= E E ((^P^l)' E E K.-K.AQpJ)- 

(19) 



Pl=0p2=0 
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When I = 0, Rq gives the joint detection probability 
that there is no hologram in both the signal and idler 
beam. And for the case I ^ 0, Ri represents the joint 
detection probability with one AZ = —I hologram in the 
signal beam and one Al = I hologram in the idler beam. 
Obviously, Ri = R-i. 

Substitute the values of G, P and Q calculated in the 
above sections to Eq. (19), we can get the relative joint 
detection probability of the three cases / = 0, 1, 2 as 



i?o : Pi : P2 = 1 : 0.346 : 0.101. 



(20) 



If we also assume that pi = p2 = for the SPDC pro- 
cess as the recent papers [22,24,25,20], the joint detection 
probability can thus be written as 

Ri = E E (p^;p'Plp'Qp'QpO^- (21) 
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And the relation for Z = 0, 1, 2 becomes: 
Ro:Ri:R2 = l: 0.311 : 0.079. 



(22) 



The difference between Eq. (20) and (22) is caused by the 

additional assumption that Lagucrrc index pi = P2 = 0. 
But this assumption has not been proven either in theo- 
retical discussions [16,20] or in experiments [22,24]. Prom 
the above calculations, we can see that this assumption 
will cause non-trivial influence to the relative joint de- 
tection probability. Our results put forward a feasible 
method to verify the assumption. 

If we rule out the influence of the holograms and mono- 
mode fibres, and make the assumption of pi = p2 = 0, the 
relationship for the relative joint detection probability of 
the cases Z = 0, 1, 2 becomes 



i?o : i?i : i?2 = 1 : 0.444 : 0.198, 



(23) 



Compare Eq. (22) with Eq. (23), we can see that the 
experimental elements can apparently influence the joint 
detection probability. 

In the experiment by Vaziri and co-workers [24], they 
found that the state of the 0AM entangled photons from 
SPDC was given by 

V' = 0.65 jO, 0) + 0.60 11, -1) + 0.47 1-1, 1) . (24) 

Prom this equation we can get the relationship for the rel- 
ative joint detection probabilities of the cases Z = 0, 1, — 1 
as 



i?o : i?i : = 1 : 0.852 : 0.523. 



(25) 



Including the influence of computer generated holograms 
and mono-mode fibres, and loosening the assumption of 
Pi = P2 = 0, we expect this relation be 



Rq : Ri : R-i = 1 : 0.346 : 0.346. 



(26) 



The reason for the difference between Eq. (25) and (26) 
might be as follows: in experiment, the diffraction ef- 
ficiency of the hologram can not be 100% . Generally, 
different holograms have different diffraction efficiencies. 
The waists of the input beam and the output beam of 
holograms will also affect the final experiment detection 
probabilities. In addition, the fibre diameter and MPD 
of the practical mono-mode fibre will also affect the de- 
tection efficiencies of different modes light. Evidently, 
the detection efficiency of avalanche detectors has little 
effect on the relative joint detection probabilities. Thus 
in the practical experiment, the P and Q values have to 
be adjusted according to the particular conditions. 



VI. CONCLUSION 

In conclusion, we have calculated the probability am- 
plitudes of different LG modes of the down converted 



photons. Our results show that the relative amplitude 
decreases almost exponentially with growing of 0AM. We 
also discussed the impact of the previous assumption for 
p on the joint detection probability. In addition, we ana- 
lyzed the influences of the experiment elements. We con- 
cluded that only a few dimensions can be explored from 
the infinite 0AM modes of the down-converted photon 
pairs. The experiment verification of the present the- 
ory is straightforward and is currently in progress in our 
laboratory. 
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Table 1. The relative probability amplitude Cp'~p^ of 

the down converted photons from SPDC. We let Cqq be 
unity. 

Fig. 1. The relative probability amplitude C for 1 = 
and pi,p2 = 0, 1,2,3,4. We let Cq q be unity. 

Fig. 2. The relative probability amplitude C for pi = 
P2 = and ? 0, 1, 2, 3, 4. We let' C^'^ be unity. 

Table 2. The relative weight of different modes after 
the computer generated hologram with M = 1. 

Table 3. The relative weight of different modes after 
the computer generated hologram with Al = 2. 
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Fig.1. 
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